One-photon double ionization of helium: a heuristic formula for the cross section 
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Without a formal derivation, we propose a formula for the total and single-differential cross in the problem 
of one-photon double ionization of an atom. The formula is benchmarked against accurate experimental data 
for the total cross section of helium. Furthermore, a direct comparison with ab initio calculations for the double 
ionization of Li + suggests that the framework is valid for the entire helium isoelectronic sequence. To this end, 
we introduce a formula for the double ionization of lithium, as well as for the triple ionization of lithium and 
beryllium. 
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Double photoionization of helium by a single photon has 
been studied for a period of more than 40 years since the pio- 
neering work of Byron and Joachain [1], who pointed out the 
importance of electron correlation in the process. Such cor- 
related processes pose many challenges, and it is only during 
the last 15 years or so that quantitative agreement between ex- 
periment H and theory IBMlotl for the total cross section was 
obtained, and a rather complete understanding of the breakup 
process in helium has emerged lfl7l - l25ll . A simple analyti- 
cal formula for the shape of single-photon multiple ionization 
cross sections was proposed by Pattard 12611 . This shape func- 
tion contains two parameters, the position and height of the 
cross section maximum. 

In this work, we propose a formula for the single- 
differential cross section in the process of one-photon double 
ionization of an atom. The formula contains a scaling factor 
that determines the height of the cross section maximum. Pro- 
vided the value of this parameter is set to one, it is shown that 
the formula yields cross sections that are in agreement with 
both theoretical and experimental double photoionization data 
for He, Li + and Li. Furthermore, it is demonstrated how the 
framework can be generalized to account for triple photoion- 
ization of lithium and beryllium. 

The present work is partly motivated by the idea behind 
a recent model for direct (nonsequential) two-photon double 
ionization of helium [27], i.e., that the explicit form of the 
electron-electron interaction is not of crucial importance in or- 
der to obtain a qualitative description of the electrons' route to 
the double continuum, as far as the total and single-differential 
cross sections are concerned. It is here argued that the corre- 
sponding one-photon double ionization event is similarly dic- 
tated by the electrons' electric dipole couplings to their re- 
spective single-particle continua, rather than the Coulombic 
interaction between the electrons. As such, the assumption is 
that the electron-electron interaction merely plays the role of 
distributing the excess energy between the ejected electrons in 
the excitation process, assuring that the total energy of the sys- 
tem is conserved. Keeping in mind that the electrons are emit- 
ted more or less simultaneously (in coincidence) in the dou- 
ble ionization process, we simply assume that there is essen- 
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tially no time for the electrons to explore the explicit geomet- 
ric form of the repulsive potential at the instant of ionization, 
which again suggests that the electron-electron interaction can 
be handled in a simple and approximate way. The simplest 
possible approximate model interaction that allows for double 
ionization by photon impact (to first order in perturbation the- 
ory), controls the energy given to each electron, and which is 
symmetric with respect to both electrons and dipole-like for 
each electron independently, can be written on the following 
heuristic form 

H int ocE(t)ziz 2 . (1) 

Here E(t) is the electric field modeling the laser pulse, which 
for simplicity is assumed to be z-polarized, and z\ and zi are 
the z-coordinates of electron 1 and 2, respectively. Although 
we give no formal proof of the assertion Eq. £[), we will nev- 
ertheless make use of it in the following to derive an explicit 
formula for the single-differential cross section in the process 
of double photoionization of helium and compare it with ac- 
curate experimental data. 

However, before proceeding, we would like to emphasize 
that the model interaction in Eq. (Q} is heuristic in nature and 
should be used with caution. For example, it fails completely 
in describing the evolution of the system in the time after the 
electrons have been emitted into the continuum, and as such 
it would in general yield incorrect angular distributions. On 
the other hand, by construction the interaction allows for the 
possibility that the electrons can absorb the photon as a unified 
system, concordant with the model of F0rre et al. If27ll . 

In the next step of approximation, both electrons are consid- 
ered to be independent particles and the ground (initial) state 
wave function of the helium atom is simply approximated by 
the product ansatz 

«P / (ri,r 2 ) = ^(r 1 )^(r 2 ), (2) 

where \j/ ls refers to the ground state of the He + -ion. It should 
be noted that the overlap between this approximate wave func- 
tion and the real ground state wave function of helium is more 
than 90%. As such, the assumption is that most of the essential 
features relevant for the double ionization process in helium 
are captured in the simplified wave function. This is a crucial 
point in the model presented here. Likewise, the final state 
wave function is approximated by a symmetrized product of 
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two He + continuum states, 
1 



*/(ri,r 2 ): 
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(3) 

with N e = 2 being the number of electrons involved in the 
ionization process. 

Applying lowest order perturbation theory to the resulting 
system, with the interaction defined in Eq. ((TJ, and the initial 
and final states defined in Eqs. (O and (01, the resulting single- 
differential cross section for double photoionization of helium 
takes the form 

^a>\(Y El \z\Wi s )\ 2 \(YE 2 \z\Wu)\ 2 , (4) 



N, 



with 



E\ +E2 = ftco — E b , 



(5) 



TiCO being the photon energy and E b = 79 eV the binding en- 
ergy of helium. As such, the total binding energy of the sys- 
tem is not considered a free parameter in the present work. 
The coupling elements in Eq. are related to the well known 
one-photon (one-electron) photoionization cross section of 
He+ JUl via the relation 129] 



a He+ oc (E-E 1s )\(We\z\Yi s 



IS 



where E\ s is the energy of the He + ground state, and <J H ^ 
the photoionization cross section of He + . 

Combining Eqs. and ©, we propose the following for- 
mula for the single-differential cross section in the process of 
one-photon double ionization of helium, 



do 
dE~i 



C 



hco <7 H , {Ei - Eis) o He + {E2 - E 



Eu 



(7) 



where Ei +E2 = fico — £),, C is a yet unknown dimensionless 
constant, and «o is the Bohr radius. Applying Eq. (0, the 
single-electron photoionization cross section of He + is multi- 
plied by two to account for the statistical weight of having two 
identical electrons in the Is orbital initially. Furthermore, the 
presence of the energy factors in the denominators is related 
to the fact that the photon is absorbed simultaneously by both 
electrons via a nonresonant transition to the final state of each 
electron 12711 . 



The one-photon double ionization process of helium has 
been investigated in length in both theoretical and experimen- 
tal studies, resulting in close quantitative agreement in the to- 
tal cross sections. In Fig. Q] we compare the total (integrated) 
cross section obtained with Eq. (|7), choosing C = 1, with the 
accurate experimental data of Samson et al. [2], who stated 
the accuracy of their results to be within ±2%. Provided we 
choose the value of the (unknown) constant C = 1, the model 
predictions is, within the experimental uncertainty, in almost 
exact agreement with the experimental data over the entire in- 
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FIG. 1 : (Color online) Double photoionization cross section of he- 
lium versus photon energy. Black line: model result Eq. with 
C = 1. Red squares: experimental results of Samson et al. |2|]. 
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FIG. 2: (Color online) Double photoionization cross section of Li + 
versus photon energy. Black line: model result. Red squares: theo- 
retical results by van der Hart and Feng [12]. Blue circles: theoretical 
results by Kheifets and Bray |7]. Green diamonds: theoretical results 
by Kleiman et al. (3(jl. 



terval of photon energies considered. The value of the param- 
eter C is now determined and will no longer be considered a 
free parameter in the rest of this work. Finally, we would like 
to note that formula (|7]i also yields single-differential cross 
sections that are in good agreement with calculated and mea- 
sured data for helium. 

Figure [2] shows the result of formula Eq. (Q when applied 
on the problem of double ionization of Li + , inserting the cor- 
responding photoionization cross sections of Li 2+ into the 
equation. Theoretical results by Kheifets and Bray |7j], van 
der Hart and Feng B12I1 and Kleiman et al. 13011 are included 
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FIG. 3: (Color online) Double photoionization cross section of 
lithium versus photon energy. Black line: model result Eq. The 
absolute one-photon photoionization cross sections of Li and Li + 
are taken from [28]. Red squares: theoretical results by Colgan et 
al. l34ll . Blue circles: experimental results by Wehlitz et al. l32ll . 
Green diamonds: experimental results by Huang et al. Bill. Black 
triangles: experimental results by Wehlitz and Juranic 13311 . 



for comparison, and the formula seems to be consistent with 
the calculated data, suggesting that it is valid for the entire 
helium isoelectronic sequence. 

In order to test the validity of the theoretical framework fur- 
ther, we now turn to the more challenging case, namely the 
double photoionization of lithium. The assumption is that it 
is the outer weakly bound 2s electron and one of the tightly 
bound Is electrons that are emitted. This suggests the follow- 
ing formula for the single-differential cross section, 



da 

dE\ 



hco 
87 



a u+ {E x -E ls ) a u {E 2 -E 2s ) 







£1 -E. 



E 2 — E 2s 



(8) 



where E\ s = —75.6 eV and E 2s = —5.4 eV are the effective 
(single-electron) energies (the negative of the ionization po- 
tential) of the Is (inner) and 2s (outer) electrons, and o Li + and 
On are the one-photon single-ionization cross sections of Li + 
and Li, respectively [28]. Figure [3] depicts the results for the 
total cross section for double ionization of lithium, as obtained 
by integrating Eq. ([SJ^and a comparison with previously ob- 
tained experimental 0TI - I33I1 and theoretical data l34ll . It turns 
out that the formula yields results that are in good agreement 
with the experimental measurements, which is somewhat sur- 
prising given the high complexity of the problem. 

Finally, we consider the problem of triple ionization of 
lithium and beryllium. We do this to show that it is relatively 
straight forward to generalize the framework to consider mul- 
tiple ionization processes. As a matter fact, combining the 
cross section for the double ionization of Li + , as obtained in 
Fig. 12 together with the photoionization cross section of neu- 
tral lithium, G[j, according to the rule in Eq. (0, the triple 
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FIG. 4: (Color online) Upper panel: Triple photoionization cross 
section of lithium versus photon energy. Black line: model re- 
sult Eq. (|9j. Green dashed line: theoretical results by Kheifets and 
Bray 13511 . Green triangles: theoretical results by Colgan et al. 13611 . 
Red diamonds, blue squares and black circles: experimental results 
by Wehlitz et al. I37H39I1 . Lower panel: Triple photoionization cross 
section of beryllium versus photon energy. Black line: model result 
Eq. U0\ . Red squares: theoretical results by Kheifets and Bray [35]. 
Blue diamonds: theoretical results by Colgan et al. l36ll . 



photoionization cross section of lithium is simply given by 

- D .{E l2 -E u+ ) a Li (E 3 -E 2s ) 



da 

dEi 



hco 



J Li 



2\a\ 



En-E A 



E-i - E 2s 



(9) 



Here o® + denotes the double photoionization cross section of 
Li + (Fig.[2]l, E u + = —198 eV is the total energy of the two 
bound Is (inner) electrons, E 2s — —5.4 eV is the energy of 
the 2s (outer) electron, £12 = E\ +E 2 , and £12 +£3 = hco + 
E Li + + E 2s is the total excess energy shared by the electrons in 
the continuum. 

The result for triple ionization of lithium is shown in Fig. [4] 
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(upper panel) and compares well with both experimental lf37l — 
[3911 and theoretical Jill H3, @3 EH data. Interestingly, the 
results are in close agreement with the predictions of the dou- 
ble shake-off model for triple photoionization, proposed by 
Kheifets and Bray 13511 . in particular for the higher photon en- 
ergies. 

The triple photoionization cross section of beryllium is cor- 
respondingly given by 

da = hco a&lEn-Ene) g gg 2+( £ 3 -E Bt s+) 
dE 3 2\a\ E n -E Be E 3 -E Be 2+ 

with a Be being the double photoionization cross section of 
Be (as calculated by Eq. (|7]l inserting the photoionization 
cross section of Be + ), cy Be i+ the one-photon single-ionization 
cross section of Be 2+ , E Be = —27.5 eV the total energy of 
the two (outer) bound 2s electrons, and E Be i+ = —153.9 eV 
the effective (single-electron) energy (the negative of the ion- 
ization potential) of the (active) Is electron. Furthermore, 
E12 =E\+ E2, and E\2 + £3 = fi(0 + E Be + E Be 2+ is the total 
excess energy shared by the three electrons in the continuum. 

The lower panel in Fig. [4] depicts the result for beryllium 
together with the theoretical calculations by Colgan et al. IT36T1 
and Kheifets and Bray 13511 . Quite interestingly, the predic- 
tion of the formula is again in favor of the result of the double 
shake-off model by Kheifets and Bray [35], as far as the to- 
tal ionization yield is concerned, but further theoretical and 
experimental investigations are required in order to settle the 
problem definitely. 

In conclusion, without a formal derivation, we have pro- 
posed a formula for the single-differential cross section in the 



problem of one-photon double ionization of He, Li + and Li. 
The corresponding function contains an unknown (dimension- 
less) constant that dictates the height of the cross section max- 
imum. The value of the constant was determined in helium by 
fitting with the experimental data of Samson et al. |2|], and 
the same value for the constant was used in the other systems. 
Provided a qualified guess for the initial state is taken, the re- 
sulting parametrization is shown to yield results in near quan- 
titative agreement with experimental and theoretical data for 
all considered cases. Finally, the problem of triple photoion- 
ization of lithium and beryllium is studied. It is demonstrated 
that agreement with experimental and theoretical results can 
be obtained. Furthermore, although not shown here, our re- 
sults are consistent with the general shape function proposed 
by Pattard [26] and the results of the half-collision model by 
Pattard and Burgdorfer B25l 14211 . 

As a final remark, we would like to add that it is rela- 
tively straight forward to generalize the framework to account 
for multiple ionization processes involving more than two or 
three electrons, like e.g. the process of quadruple ionization of 
beryllium, which is a problem that is difficult to pursue within 
the framework of more rigorous ab initio methods. 
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